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In a closed economic system, money is conserved. Thus, by analogy with energy, the equilibrium 
probability distribution of money must follow the exponential Boltzmann-Gibbs law characterized by 
an effective temperature equal to the average amount of money per economic agent. We demonstrate 
how the Boltzmann-Gibbs distribution emerges in computer simulations of economic models. Then 
we consider a thermal machine, in which the difference of temperatures allows one to extract a 
monetary profit. We also discuss the role of debt, and models with broken time-reversal symmetry 
for which the Boltzmann-Gibbs law does not hold. The instantaneous distribution of money among 
the agents of a system should not be confused with the distribution of wealth. The latter also 
includes material wealth, which is not conserved, and thus may have a different (e. g. power-law) 
distribution. 



PACS numbers: 87.23. Ge, 05.90.+m, 89.90. +n, 02.50.- 



I. INTRODUCTION 

The application of statistical physics methods to eco- 
nomics promises fresh insights into problems tradition- 
ally not associated with physics (see, for example, the 
recent review and book Both statistical mechanics 
and economics study big ensembles: collections of atoms 
or economic agents, respectively. The fundamental law of 
equilibrium statistical mechanics is the Boltzmann-Gibbs 
law, which states that the probability distribution of en- 
ergy e is P{e) = Ce~ e l T , where T is the temperature, 
and C is a normalizing constant Q. The main ingredi- 
ent that is essential for the textbook derivation of the 
Boltzmann-Gibbs law Q is the conservation of energy 
||. Thus, one may generalize that any conserved quan- 
tity in a big statistical system should have an exponential 
probability distribution in equilibrium. 

An example of such an unconventional Boltzmann- 
Gibbs law is the probability distribution of forces expe- 
rienced by the beads in a cylinder pressed with an ex- 
ternal force Q. Because the system is at rest, the total 
force along the cylinder axis experienced by each layer of 
granules is constant and is randomly distributed among 
the individual beads. Thus the conditions are satisfied 
for the applicability of the Boltzmann-Gibbs law to the 
force, rather than energy, and it was indeed found exper- 
imentally ||]. 

We claim that, in a closed economic system, the to- 
tal amount of money is conserved. Thus the equilib- 
rium probability distribution of money P(m) should fol- 
low the Boltzmann-Gibbs law P(m) — Ce~ m l T . Here 
to is money, and T is an effective temperature equal to 
the average amount of money per economic agent. The 
conservation law of money j^] reflects their fundamen- 
tal property that, unlike material wealth, money (more 
precisely the fiat, "paper" money) is not allowed to be 
manufactured by regular economic agents, but can only 
be transferred between agents. Our approach here is very 



similar to that of Ispolatov et al. However, they con- 
sidered only models with broken time-reversal symmetry, 
for which the Boltzmann-Gibbs law typically does not 
hold. The role of time-reversal symmetry and deviations 
from the Boltzmann-Gibbs law are discussed in detail in 



Sec. VII 



It is tempting to identify the money distribution P(m) 
with the distribution of wealth ||. However, money 
is only one part of wealth, the other part being ma- 
terial wealth. Material products have no conservation 
law: They can be manufactured, destroyed, consumed, 
etc. Moreover, the monetary value of a material product 
(the price) is not constant. The same applies to stocks, 
which economics textbooks explicitly exclude from the 
definition of money So, in general, we do not ex- 
pect the Boltzmann-Gibbs law for the distribution of 
wealth. Some authors believe that wealth is distributed 
according to a power law (Pareto-Zipf), which originates 
from a multiplicative random process g). Such a pro- 
cess may reflect, among other things, the fluctuations of 
prices needed to evaluate the monetary value of material 
wealth. 



II. BOLTZMANN-GIBBS DISTRIBUTION 

Let us consider a system of many economic agents 
N 3> 1, which may be individuals or corporations. In 
this paper, we only consider the case where their number 
is constant. Each agent i has some money and may 
exchange it with other agents. It is implied that money is 
used for some economic activity, such as buying or selling 
material products; however, we are not interested in that 
aspect. As in Ref. ||, for us the only result of interaction 
between agents i and j is that some money Am changes 



hands: [mi,rrij] 



j] = [rm — Am, rrij + Am] . No- 



tice that the total amount of money is conserved in each 
transaction: to; + m,- = to' + to' . This local conserva- 
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tion law of money || is analogous to the conservation 
of energy in collisions between atoms. We assume that 
the economic system is closed, i. e. there is no external 
flux of money, thus the total amount of money M in the 
system is conserved. Also, in the first part of the paper, 
we do not permit any debt, so each agent's money must 
be non-negative: rrii > 0. A similar condition applies to 
the kinetic energy of atoms: £j > 0. 

Let us introduce the probability distribution function 
of money P(m), which is defined so that the number of 
agents with money between m and m + dm is equal to 
NP(m) dm. We are interested in the stationary distribu- 
tion P(m) corresponding to the state of thermodynamic 
equilibrium. In this state, an individual agent's money 
rrii strongly fluctuates, but the overall probability distri- 
bution P(m) does not change. 

The equilibrium distribution function P(m) can be de- 
rived in the same manner as the equilibrium distribution 
function of energy P(e) in physics 0]. Let us divide the 
system into two subsystems 1 and 2. Taking into account 
that money is conserved and additive: m — m\ + m,2, 
whereas the probability is multiplicative: P = P1P2, 
we conclude that P{m\ + m 2 ) = P(mi)P(m,2). The 
solution of this equation is P(m) = Ce~ m l T \ thus the 
equilibrium probability distribution of money has the 
Boltzmann-Gibbs form. From the normalization condi- 
tions J °° P(m) dm — 1 and / °° mP(m) dm — M/N, we 
find that C = 1/T and T = M/N. Thus, the effective 
temperature T is the average amount of money per agent. 

The Boltzmann-Gibbs distribution can be also ob- 
tained by maximizing the entropy of money distribu- 
tion S = — J °° dm P(m) In P(m) under the constraint of 
money conservation [^j. Following original Boltzmann's 
argument, let us divide the money axis < m < 00 
into small bins of size dm and number the bins con- 
secutively with the index b = 1,2,... Let us denote 
the number of agents in a bin b as Nb, the total num- 
ber being N — Ylb=i ^b- The agents in the bin b have 
money nib, an d the total money is M — Y^=i m bNb- 
The probability of realization of a certain set of occu- 
pation numbers {Nb} is proportional to the number of 
ways iV agents can be distributed among the bins pre- 
serving the set {N b }. This number is Nl/NilN 2 \ ■ ■ . The 
logarithm of probability is entropy In AH — Ylb=i hiJV(,!. 
When the numbers A'f, are big and Stirling's formula 
In AT! w N\nN applies, the entropy per agent is S = 
(NlnN - EZi N b^N b )/N = -E^iAlnP 6 , where 
Pb = Nb/N is the probability that an agent has money 
TOf,. Using the method of Lagrange multipliers to maxi- 
mize the entropy 5* with respect to the occupation num- 
bers {Nb} with the constraints on the total money M and 
the total number of agents N generates the Boltzmann- 
Gibbs distribution for P(m) M. 



N=500, M=5*10 5 , time=4*10 5 . 
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FIG. 1. Histogram and points: stationary probability dis- 
tribution of money P(ra). Solid curves: fits to the Boltz- 
mann-Gibbs law P(m) oc exp(— m/T). Vertical lines: the 
initial distribution of money. 

III. COMPUTER SIMULATIONS 

To check that these general arguments indeed work, 
we performed several computer simulations. Initially, all 
agents are given the same amount of money: P(m) — 
S(m — M/N), which is shown in Fig. [l] as the double 
vertical line. One pair of agents at a time is chosen ran- 
domly, then one of the agents is randomly picked to be 
the "winner" (the other agent becomes the "loser" ) , and 
the amount Am > is transferred from the loser to the 
winner. If the loser does not have enough money to pay 
{mi < Am), then the transaction does not take place, 
and we proceed to another pair of agents. Thus, the 
agents are not permitted to have negative money. This 
boundary condition is crucial in establishing the station- 
ary distribution. As the agents exchange money, the 
initial delta-function distribution first spread symmetri- 
cally. Then, the probability density starts to pile up at 
the impenetrable boundary m = 0. The distribution be- 
comes asymmetric (skewed) and ultimately reaches the 
stationary exponential shape shown in Fig. [l]. We used 
several trading rules in the simulations: the exchange 
of a small constant amount Am = 1, the exchange of 
a random fraction < v < 1 of the average money of 
the pair: Am = i/(m, + m,j)/2, and the exchange of a 
random fraction v of the average money in the system: 
Am = v M/N. Figures in the paper mostly show sim- 
ulations for the third rule; however, the final stationary 
distribution was found to be the same for all rules. 

In the process of evolution, the entropy S increases 
in time and saturates at the maximal value for the 
Boltzmann-Gibbs distribution. This is illustrated by the 
top curve in Fig. |^ computed for the third rule of ex- 
change. The bottom curve in Fig. p] shows the time evo- 
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FIG. 2. Time evolution of entropy. Top curve: for the 
exchange of a random fraction v of the average money in the 
system: Am = vM/N. Bottom curve: for the exchange 
of a small constant amount Am = 1. The time scale for 
the bottom curve is 500 times greater than indicated, so it 
actually ends at the time 10 6 . 

lution of entropy for the first rule of exchange. The time 
scale for this curve is 500 times greater than for the top 
curve, so the bottom curve actually ends at the time 10 6 . 
The plot shows that, for the first rule of exchange, mixing 
is much slower than for the third one. Nevertheless, even 
for the first rule, the system also eventually reaches the 
Boltzmann-Gibbs state of maximal entropy, albeit over a 
time much longer than shown in Fig. ||. 

One might argue that the pairwise exchange of money 
may correspond to a medieval market, but not to a mod- 
ern economy. In order to make the model somewhat more 
realistic, we introduce firms. One agent at a time be- 
comes a "firm" . The firm borrows capital K from an- 
other agent and returns it with an interest rK, hires L 
agents and pays them wages W, manufactures Q items of 
a product and sell it to Q agents at a price R. All of these 
agents are randomly selected. The firm receives the profit 
F = RQ — LW — rK. The net result is a many-body ex- 
change of money that still satisfies the conservation law. 

Parameters of the model are selected following the pro- 
cedure described in economics textbooks. The aggre- 
gate demand-supply curve for the product is taken to 
be R{Q) — V/Q n , where Q is the quantity people would 
buy at a price i?, and T) = 0.5 and V = 100 are constants. 
The production function of the firm has the conventional 
Cobb-Douglas form: Q(L, K) = L^K 1 '^, where (3 = 0.8 
is a constant. In our simulation, we set W = 10. By 
maximizing firm's profit F with respect to K and L, we 
find the values of the other parameters: L = 20, Q = 10, 
R = 32, and F = 68. 

However, the actual values of the parameters do not 
matter. Our computer simulations show that the sta- 
tionary probability distribution of money in this model 



always has the universal Boltzmann-Gibbs form indepen- 
dent of the model parameters. 

IV. THERMAL MACHINE 

As explained in Introduction, the money distribution 
P(m) should not be confused with the distribution of 
wealth. We believe that P(m) should be interpreted as 
the instantaneous distribution of purchasing power in the 
system. Indeed, to make a purchase, one needs money. 
Material wealth normally is not used directly for a pur- 
chase. It needs to be sold first to be converted into 
money. 

Let us consider an outside monopolistic vendor sell- 
ing a product (say, cars) to the system of agents at a 
price p. Suppose that a certain small fraction / of the 
agents needs to buy the product at a given time, and 
each agent who has enough money to afford the price 
will buy one item. The fraction / is assumed to be suffi- 
ciently small, so that the purchase does not perturb the 
whole system significantly. At the same time, the ab- 
solute number of agents in this group is assumed to be 
big enough to make the group statistically representative 
and characterized by the Boltzmann-Gibbs distribution 
of money. The agents in this group continue to exchange 
money with the rest of the system, which acts as a ther- 
mal bath. The demand for the product is constantly 
renewed, because products purchased in the past expire 
after a certain time. In this situation, the vendor can 
sell the product persistently, thus creating a small steady 
leakage of money from the system to the vendor. 

What price p would maximize the vendor's income? 
To answer this question, it is convenient to introduce 
the cumulative distribution of purchasing power M{m) = 
N P(m') dm' = Ne~ m ^ T , which gives the number of 
agents whose money is greater than m. The vendor's 
income is fpj\f(p). It is maximal when p = T, i. e. the 
optimal price is equal to the temperature of the system. 
This conclusion also follows from the simple dimensional 
argument that temperature is the only money scale in the 
problem. At the price p — T that maximizes the vendor's 
income, only the fraction N{T)/N — e" 1 — 0.37 of the 
agents can afford to buy the product. 

Now let us consider two disconnected economic sys- 
tems, one with the temperature T\ and another with T^: 
T\ > Ti. A vendor can buy a product in the latter sys- 
tem at its equilibrium price Ti and sell it in the former 
system at the price Ti, thus extracting the speculative 
profit Ti — T-2, as in a thermal machine. This example 
suggests that speculative profit is possible only when the 
system as a whole is out of equilibrium. As money is 
transferred from the high- to the low-temperature sys- 
tem, their temperatures become closer and eventually 
equal. After that, no speculative profit is possible, which 
would correspond to the "thermal death" of the economy. 
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N=500, M=5*10 5 , time=4*10 5 . 
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FIG. 3. Histograms: stationary distributions of money 
with and without debt. Solid curves: fits to the Boltz- 
mann-Gibbs laws with temperatures T = 1800 and T — 1 000. 

This example brings to mind economic relations between 
developed and developing countries, with manufacturing 
in the poor (low-temperature) countries for export to the 
rich (high-temperature) ones. 

V. MODELS WITH DEBT 

Now let us discuss what happens if the agents are per- 
mitted to go into debt. Debt can be viewed as negative 
money. Now when a loser does not have enough money to 
pay, he can borrow the required amount from a reservoir, 
and his balance becomes negative. The conservation law 
is not violated: The sum of the winner's positive money 
and loser's negative money remains constant. When an 
agent with a negative balance receives money as a win- 
ner, she uses this money to repay the debt until her bal- 
ance becomes positive. We assume for simplicity that the 
reservoir charges no interest for the lent money. However, 
because it is not sensible to permit unlimited debt, we put 
a limit md on the maximal debt of an agent: m, > —md- 
This new boundary condition P(m < —m e i) = replaces 
the old boundary condition P(m < 0) = 0. The result of 
a computer simulation with md = 800 is shown in Fig. 
[I together with the curve for md = 0. P(m) is again 
given by the Boltzmann-Gibbs law, but now with the 
higher temperature T = M/N + md, because the nor- 
malization conditions need to be maintained including 
the population with negative money: f^° m P(m) dm = 1 
and mP(m) dm = M/N. The higher temperature 
makes the money distribution broader, which means that 
debt increases inequality between agents |l. 

Imposing a sharp cutoff at md may be not quite real- 
istic. In practice, the cutoff may be extended over some 



range depending on the exact bankruptcy rules. Over 
this range, the Boltzmann-Gibbs distribution would be 
smeared out. So we expect to see the Boltzmann-Gibbs 
law only sufficiently far from the cutoff region. Similarly, 
in experiment [Q, some deviations from the exponential 
law were observed near the lower boundary of the distri- 
bution. Also, at the high end of the distributions, the 
number of events becomes small and statistics poor, so 
the Boltzmann-Gibbs law loses applicability. Thus, we 
expect the Boltzmann-Gibbs law to hold only for the in- 
termediate range of money not too close either to the 
lower boundary or to the very high end. However, this 
range is the most relevant, because it covers the great 
majority of population. 

Lending creates equal amounts of positive (asset) and 
negative (liability) money When economics text- 

books describe how "banks create money" or "debt cre- 
ates money" 0, they do not count the negative liabil- 
ities as money, and thus their money is not conserved. 
In our operational definition of money, we include all 
financial instruments with fixed denomination, such as 
currency, IOUs, and bonds, but not material wealth or 
stocks, and we count both assets and liabilities. With 
this definition, money is conserved, and we expect to see 
the Boltzmann-Gibbs distribution in equilibrium. Unfor- 
tunately, because this definition differs from economists' 
definitions of money (Ml, M2, M3, etc. @]), it is not 
easy to find the appropriate statistics. Of course, money 
can be also emitted by a central bank or government. 
This is analogous to an external influx of energy into a 
physical system. However, if this process is sufficiently 
slow, the economic system may be able to maintain quasi- 
equilibrium, characterized by a slowly changing temper- 
ature. 

We performed a simulation of a model with one bank 
and many agents. The agents keep their money in ac- 
counts on which the bank pays interest. The agents may 
borrow money from the bank, for which they must pay 
interest in monthly installments. If they cannot make 
the required payments, they may be declared bankrupt, 
which relieves them from the debt, but the liability is 
transferred to the bank. In this way, the conservation 
of money is maintained. The model is too elaborate to 
describe it in full detail here. We found that, depend- 
ing on the parameters of the model, either the agents 
constantly lose money to the bank, which steadily re- 
duces the agents' temperature, or the bank constantly 
loses money, which drives down its own negative balance 
and steadily increases the agents' temperature. 

VI. BOLTZMANN EQUATION 

The Boltzmann-Gibbs distribution can be also derived 
from the Boltzmann equation [n0|, which describes the 
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FIG. 4. Histogram: stationary probability distribution of 
money in the multiplicative random exchange model studied 
in Ref. [HI. Solid curve: the Boltzmann-Gibbs law. 



N=500, M=5*10 5 , tax=40%. 
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FIG. 5. Histogram: stationary probability distribution of 
money in the model with taxes and subsidies. Solid curve: 
the Boltzmann-Gibbs law. 



time evolution of the distribution function P(m) due to 
pairwise interactions: 



dP(m) 
dt 



{—w\ +A]P(m)P(m') (1) 



+W[ m -A,m'+A]^[ m , m ']P(m - A)P(m' + A)} dm' dA. 

Here u>[ m . m /]^[ m _A,m'+A] is the rate of transferring 
money A from an agent with money to to an agent 
with money to'. If a model has time-reversal symme- 
try, then the transition rate of a direct process is the 
same as the transition rate of the reversed process, thus 
the w-factors in the first and second lines of Eq. ([j]) are 
equal. In this case, the Boltzmann-Gibbs distribution 
P(m) = C exp(— to/T) nullifies the right-hand side of Eq. 
(|l)); thus this distribution is stationary: dP(m)/dt = 



VII. NON-BOLTZMANN-GIBBS 
DISTRIBUTIONS 

However, if time-reversal symmetry is broken, the two 
transition rates w in Eq. (p]) may be different, and the 
system may have a non-Boltzmann-Gibbs stationary dis- 
tribution or no stationary distribution at all. Exam- 
ples of such kind were studied in Ref. 0. One model 
was called the multiplicative random exchange. In this 
model, a randomly selected loser i loses a fixed frac- 
tion a of his money to a randomly selected winner 
j: [mi^rnj] — > [(1 — a)mi , m.j + ami]. If we try to 
reverse this process and appoint the winner j to be- 
come a loser, the system does not return to the origi- 
nal configuration [to^to.,]: [(1 — a)mi , rrij + ami] — > 
[(1 — a)rrii + aijiij + ami) , (1 — a)(nij + ami)]. Except 
for a — 1/2, the exponential distribution function is not a 



stationary solution of the Boltzmann equation derived for 
this model in Ref. Instead, the stationary distribu- 
tion has the shape shown in Fig. |]for a = 1/3, which we 
reproduced in our numerical simulations. It still has an 
exponential tail end at the high end, but drops to zero at 
the low end for a < 1/2. Another example of similar kind 
was studied in Ref. pl|, which appeared after the first 



version of our paper was posted as cond-mat/0001432 
on January 30, 2000. In that model, the agents save 
a fraction A of their money and exchange a random 
fraction e of their total remaining money: [to^TOj] — > 
[Ami + e(l — A)(toj + mj) , Amj + (1 — e)(l — \){mi+m.j)]. 
This exchange also does not return to the original config- 
uration after being reversed. The stationary probability 
distribution was found in Ref. Jll| to be nonexponential 
for A / with a shape qualitatively similar to the one 
shown in Fig. ^. 

Another interesting example of a non-Boltzmann- 
Gibbs distribution occurs in a model with taxes and sub- 
sidies. Suppose a special agent ("government") collects 
a fraction ("tax") of every transaction in the system. 
The collected money is then equally divided between all 
agents of the system, so that each agent receives the sub- 
sidy 8m with the frequency 1/t s . Assuming that 8m is 
small and approximating the collision integral with a re- 
laxation time T r flO[| , we obtain the following Boltzmann 
equation 



dP(m) 8m dP{m) _ P(m) - P(m) 



dm 



(2) 



where P(m) is the equilibrium Boltzmann-Gibbs func- 
tion. The second term in the left-hand side of Eq. (|^) 
is analogous to the force applied to electrons in a metal 
by an external electric field p0| . The approximate sta- 
tionary solution of Eq. (0) is the displaced Boltzmann- 
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Gibbs one: P(m) = P(m — (r r /r s ) 5m). The displace- 
ment of the equilibrium distribution P(m) by (r r /r s ) 5m 
would leave an empty gap near m = 0. This gap is 
filled by interpolating between zero population at m = 
and the displaced distribution. The curve obtained in a 
computer simulation of this model (Fig. ||) qualitatively 
agrees with this expectation. The low-money popula- 
tion is suppressed, because the government, acting as an 
external force, "pumps out" that population and pushes 
the system out of thermodynamic equilibrium. We found 
that the entropy of the stationary state in the model with 
taxes and subsidies is few percents lower than without. 

These examples show that the Boltzmann- Gibbs dis- 
tribution is not fully universal, meaning that it does 
not hold for just any model of exchange that conserves 
money. Nevertheless, it is universal in a limited sense: 
For a broad class of models that have time-reversal sym- 
metry, the stationary distribution is exponential and does 
not depend on the details of a model. Conversely, when 
time-reversal symmetry is broken, the distribution may 
depend on model details. The difference between these 
two classes of models may be rather subtle. For exam- 
ple, let us change the multiplicative random exchange 
from a fixed fraction of loser's money to a fixed fraction 
of the total money of winner and loser. This modifica- 
tion retains the multiplicative idea that the amount ex- 
changed is proportional to the amount involved, but re- 
stores time-reversal symmetry and the Boltzmann-Gibbs 
distribution. In the model with Am = 1 discussed in 
the next Section, the difference between time-reversible 
and time-irreversible formulations amounts to the differ- 
ence between impenetrable and absorbing boundary con- 
ditions at m — 0. Unlike in physics, in economy there is 
no fundamental requirement that interactions have time- 
reversal symmetry. However, in the absence of detailed 
knowledge of real microscopic dynamics of economic ex- 
change, the semiuniversal Boltzmann-Gibbs distribution 
appears to be a natural starting point. 

Moreover, deviations from the Boltzmann-Gibbs law 
may occur only if the transition rates w in Eq. (^) ex- 
plicitly depend on the agents money m or m' in an asym- 
metric manner. In another simulation, we randomly pre- 
selected winners and losers for every pair of agents 
In this case, money flows along directed links between 
the agents: i — > j — > k, and time-reversal symmetry is 
strongly broken. This model is closer to the real econ- 
omy, in which, for example, one typically receives money 
from an employer and pays it to a grocer, but rarely 
the reverse. Nevertheless, we still found the Boltzmann- 
Gibbs distribution of money in this model, because the 
transition rates w do not explicitly depend on m and m'. 



VIII. NONLINEAR BOLTZMANN EQUATION 
VS. LINEAR MASTER EQUATION 

For the model where agents randomly exchange the 
constant amount Am = 1, the Boltzmann equation is: 



dP„ 
dt 



Pm+1 ^ Pn Pm — 1 ^ P71 



71=0 



-p Vp -p Vp 



(3) 



n=0 



— (Pm+1 + Pm-1 — 2P m ) + Po{P m — Pm-l), (4) 

where P m = P{m) and we have used X)m=o ^ >m = 1- ^he 
first, diffusion term in Eq. (^) is responsible for broaden- 
ing of the initial delta-function distribution. The second 
term, proportional to Po, is essential for the Boltzmann- 
Gibbs distribution P m = e~ m / T (l - e~ 1/T ) to be a sta- 
tionary solution of Eq. (^). In a similar model studied in 
Ref . , the second term was omitted on the assumption 
that agents who lost all money are eliminated: Po = 0. 
In that case, the total number of agents is not conserved, 
and the system never reaches any stationary distribution. 
Time-reversal symmetry is violated, since transitions into 
the state m = are permitted, but not out of this state. 

If we treat Pq as a constant, Eq. (Q) looks like a linear 
Fokker-Planck equation jl0| for P TO , with the first term 
describing diffusion and the second term an external force 
proportional to Pq. Similar equations were studied in 
Ref. (D . Eq. (||) can be also rewritten as 



dP m 

dt 



P m+1 - (2 - P )P m + (1 - P )P m _i. (5) 



The coefficient (1 — Pq) in front of P m -i represents the 
rate of increasing money by Am — 1, and the coefficient 1 
in front of P m +i represents the rate of decreasing money 
by Am = — 1. Since Pq > 0, the former is smaller than 
the latter, which results in the stationary Boltzmann- 
Gibbs distributions P m = (I — Po) m . An equation similar 
to Eq. (||) describes a Markov chain studied for strate- 
gic market games in Ref. fl2fl . Naturally, the station- 
ary probability distribution of wealth in that model was 
found to be exponential |l2| . 

Even though Eqs. and (||) look like linear equa- 
tions, nevertheless the Boltzmann equation (Q) and (||) is 
a profoundly nonlinear equation. It contains the product 
of two probability distribution functions P in the right- 
hand side, because two agents are involved in money 
exchange. Most studies of wealth distribution ||] have 
the fundamental flaw that they use a single-particle ap- 
proach. They assume that the wealth of an agent may 
change just by itself and write a linear master equation 
for the probability distribution. Because only one parti- 
cle is considered, this approach cannot adequately incor- 
porate conservation of money. In reality, an agent can 
change money only by interacting with another agent, 
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thus the problem requires a two-particle probability dis- 
tribution function. Using Boltzmann's molecular chaos 
hypothesis, the two-particle function is factorized into 
a product of two single-particle distributions functions, 
which results in the nonlinear Boltzmann equation. Con- 
servation of money is adequately incorporated in this 
two-particle approach, and the universality of the expo- 
nential Boltzmann-Gibbs distribution is transparent. 



IX. CONCLUSIONS 

Everywhere in the paper we assumed some randomness 
in the exchange of money. The results of our paper would 
apply the best to the probability distribution of money in 
a closed community of gamblers. In more traditional eco- 
nomic studies, the agents exchange money not randomly, 
but following deterministic strategies, such as maximiza- 
tion of utility functions f^,|l3|| . The concept of equilibrium 
in these studies is similar to mechanical equilibrium in 
physics, which is achieved by minimizing energy or max- 
imizing utility. However, for big ensembles, statistical 
equilibrium is a more relevant concept. When many het- 
erogeneous agents deterministically interact and spend 
various amounts of money from very little to very big, 
the money exchange is effectively random. In the future, 
we would like to uncover the Boltzmann-Gibbs distribu- 
tion of money in a simulation of a big ensemble of eco- 
nomic agents following realistic deterministic strategies 
with money conservation taken into account. That would 
be the economics analog of molecular dynamics simula- 
tions in physics. While atoms collide following fully de- 
terministic equations of motion, their energy exchange is 
effectively random due to complexity of the system and 
results in the Boltzmann-Gibbs law. 

We do not claim that the real economy is in equilib- 
rium. (Most of the physical world around us is not in 
true equilibrium either.) Nevertheless, the concept of 
statistical equilibrium is a very useful reference point for 
studying nonequilibrium phenomena. 

The authors are grateful to M. Gubrud for helpful 
discussion and proofreading of an earlier version of the 
manuscript. 

Note added: After the paper had been submitted for 
publication, we have learned about the book by Aoki 
fL4j , who applied many ideas of statistical physics to 
economics, albeit not specifically to the distribution of 
money. 
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